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Abstract. We consider in this paper the wcll-posedness for the Cauchy prob- 
lem associated to two-dimensional dispersive systems of Boussinesq type which 
model weakly nonlinear long wave surface waves. We emphasize the case of 
the strongly dispersive ones with focus on the "KdV-KdV" system which pos- 
sesses the strongest dispersive properties and which is a vector two-dimensional 
extension of the classical KdV equation. 

1. Introduction 

1.1. General Setting. In general, nonlinear dispersive equations and systems are 
not derived from first principles. They are obtained as asymptotic models (normal 
forms) when some small parameter tends to zero, of more general systems, under 
a suitable scaling. There are supposed to describe the dynamics under suitable 
scaling conditions. We will study in this article small amplitude, long wave models 
which have the general form 

(1) U t +AU + eF{U, VU) + eLU = 0. 

Here e is a "small" parameter which takes into account the nonlinear and dis- 
persive effects, which are assumed of the same order. U — U(x, t) is a vector in 
l n+1 , n = 1, 2 and x e M", t > 0. The order zero part U t + AU is just the linear 
wave equation, while the third and fourth terms in ([T]) represents the nonlinear and 
dispersive effects. The right-hand side of JTJ should in fact be 0(e 2 ). 

A typical example is that of the abed Boussinesq systems for long, small ampli- 
tude gravity surface water waves introduced in [BJ, [5] 



(2) 



rjt + V • v + e[V • (rjv) + aV • Av - bAr) t ] = 
v 4 + V?y + e[±V|v| 2 + cVAt? - dAv t ] = 0. 

Here, r\ is the deviation of the free surface from its rest state and v is an 0(e 2 ) 
approximation of the horizontal velocity taken at a certain depth (see [B], [5]). The 
constants a, 6, c, d are modeling parameters subject to the constraint a+b+c+d = -|. 
Those three degrees of freedom arise from the choice of the height at which the 
velocity is taken and from a double use of the BBM trick. 

The small parameter e is defined by 

e = a/h~ (h/X) 2 , 
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where h denotes the mean depth of the fluid, a a typical amplitude of the wave 
and A a typical horizontal wavelength. 

Those systems are approximations of the full water wave system in the so-called 
Boussinesq regime (see [6], [5]). They degenerate into the KdV (or BBM) equation 
in the one-dimensional case, for waves traveling in one direction. They also appear 
as models for internal waves, in an appropriate regime (see [9], When surface 

tension effects are taken into account, the coefficient c should be changed into c — r 
where r > is the surface tension parameter (see [H]). Throughout the paper we 
will consider only the case of purely gravity waves, that is r = 0. 

Of course restrictions are to be imposed on a, b, c, d in order that the linear part 
of ((21) be well-posed. It was established in [7] that, when n = 1, all the linearly 
well-posed systems are locally nonlinearly well-posed. As for the two-dimensional 
case (n = 2), it has been proved in |15j that in the generic case where b > 0, d > 
one has well-posedness on time scales of order 0{-^=) for data in the Sobolev space 

H . The local well-posedness in other cases (but not in the strongly dispersive 
"KdV-KdV" case) was proved in [13] , but the question of the dependence of the 
time existence with respect to e was not addressed there. 

One of the goals of the present paper is to complete the local well-posedness 
theory for a, b, c, d systems in two dimensions, with focus on the size of the lifespan 
of the solutions with respect to e. 

An interesting fact (already noticed in [6]) from the PDE view point is that, 
though the a, b, c, d systems describe the same wave propagation phenomenon, their 
dispersive properties are quite different. We will precise the dispersion matrix in 
Section 2. After diagonalization, one is led to two-dimensional systems with strong 
dispersion (of KdV or Schrodinger type) or with weak dispersion (for instance of 
BBM type). One also obtains a system which can be viewed as a dispersive pertur- 
bation of the two-dimensional Saint- Venant (shallow-water) system, generalizing 
the one-dimensional system studied by Amick [5] and Schonbck |30j . 

On the other hand, the solutions of the Cauchy problem associated to Q should 
exist on time scales of order O(^) in order to prove that the asymptotic system is 
a good approximation of the water wave system [T]. In fact, the solutions of 
cease to be relevant as approximations of those of the original system on time scales 
larger than O(^). 

It is worth noticing that it is unlikely that all the abed systems would have 
global solutions. In fact they are shown to be Hamiltonian (and thus possess a 
formally conserved energy) only when b = d (see [6]). This can be used in the 
one-dimensional case to obtain the global well-posedness of the Cauchy problem 
for a few of the systems see [7] but this situation is exceptional and moreover does 
not apply to the two-dimensional case where the Hamiltonian does not control any 
Sobolev type norm. 

This situation is in strong contrast with one- directional models such as the 
Korteweg- de Vries (KdV) or Benjamin-Bona-Mahony (BBM) equations or with 
quasi one- directional models such as the Kadomtsev-Petviashvili (KP) equation 
where global a priori bounds are available (possibly at low degree of regularity), 
which allows to prove the global well-posedness of the Cauchy problem. 
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As far as we know the problem of solving any of the systems ([2]) on time intervals 
of order 0{\) has been completely open until very recently, except for some one- 
dimensional ones which have global solutions (under some restrictions on the initial 
data). 

In [5] Bona , Colin and Lannes obtained an order O(^) existence interval for a 
fully symmetric class of systems, in one and two dimensions. Those systems have a 
skew-adjoint linear part (a = c) and arc obtained from ((2|) via the nonlinear change 
of variables v = v(l + ^rf) which does not affect the linear part (modulo higher 
order terms in e) and which symmetrizes the nonlinear part. Neglecting the higher 
order terms in e, one obtain skew-adjoint perturbations of symmetric quasilinear 
hyperbolic systems and the classical theory of this kind of systems provides the 
O(-) existence time. This method does not use the dispersive part of the systems 
and of course does not solve the long time existence problem for the abed systems. 

In [37] an approach based on a Nash-Moser theorem is developed to prove wcll- 
posedness results for the 2D- Boussinesq systems on time intervals of order 1/e. 
In [27] the generic case, that is when 6 > 0, d > 0,a < 0,c < and the BBM/BBM 
case, that is when 6>0,<i>0,a = c = are emphasized but the approach could 
very likely be applied to other cases . This method does not use the dispersive part 
of the systems and a high regularity level in required on the initial data, with loss 
of regularity on the solution. 

The aim of the present paper is to prove the well-posedness for the most disper- 
sive of the two-dimensional Boussinesq systems ([2]) on time scales of order e -1 / 2 , 
using the dispersive properties of the systems. This allows to consider relatively 
rough initial data. We will in particular obtain uniformly bounded in e solutions on 
time scales of order T/y/e in suitable Sobolev spaces, achieving the rigorous justi- 
fication of those models on the corresponding time scales (see [1]). Note however 
that we are not able in our functional setting to reach the optimal time scales 0(-). 

The heart of the paper concerns the most dispersive Boussinesq system (the KdV- 
KdV system) where a = c = 1/6, b = d = which is an interesting two-dimensional 
extension of the Korteweg-de Vries equation for which the local well-posedness 
is not a simple matter. Our method deeply lies on dispersive estimates for the 
underlying linear problem. In particular we establish new Strichartz and maximal 
function estimates. 

For the other cases (which are less dispersive), at least when d > 0,0 the proofs 
are the extension to the two-dimensional case of those given by energy estimates 
in [7J for the one-dimensional case, keeping track of the e's but we do not pursue 
this issue here. We will focus instead on the order 2 Boussinesq systems which can 
be written as systems of nonlinear nonlocal Schrodinger type equations coupled 
by the nonlinear terms. Some of those cases have been studied in [T5] but the e 
dependence of the existence time interval was not addressed here and the difficulty 
linked to the case d = was underestimated. 

1.2. Organization of the paper. The paper is organized as follows. 
Section 2 recall briefly some useful facts on the Boussinesq systems. 



This approach gives also similar results in the one dimensional case. 
2 There are relatively few physically relevant systems of this form. 
When d = it is unclear whether the local well-posedness can be obtained by elementary 
energy methods in the 2D case without imposing the irrotationality of the velocity. 
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In Section 3 we will consider the KdV-KdV Boussinesq system, and establish the 
well-posedness in the two-dimensional case, on time intervals of order e -1 ' 2 , after 
establishing various dispersive estimates for the solutions of the linear part. 

In Section 4, we consider the well-posedness of the remaining strongly dispersive 
cases ( Schrodingcr type) to prove that there are also wcll-poscd on time intervals 
of order e -1 / 2 . There are essentially two cases, depending whether b or d vanishes. 

1.3. Notations. For any positive numbers a and b, the notation a < 6 means that 
there exists a positive constant c such that a < cb. C or c will also denote various 
positive constants independent of e. 

We denote the horizontal variables by x when n = 1 and by x = (2:1,2:2) when 
n = 2. We will also denote by • the euclidian scalar product of two vectors x = 
(x\,X2) and y = (2/1 , 2/2) 01 K 2 , which is to say x ■ y = x\y\ + X2IJ2, and by \x\ the 
euclidian norm of x = (2:1,2:2), *- e - M = \/ x \ + x \- 

We use the Fourier multiplier notation: f(D)u is defined as .F(/(.D)u)(£) = 
/(£)w(£), where T and^stand for the Fourier transform, which is defined by 

nm) = m = ^ / <Kx)e-***dx, 

for any function <fi : M 2 — >• C. 

For s £ K, wc define the Bessel and Ricsz potentials of order — s by A s = 
(1 + |_D| 2 ) S / 2 and D s = \D\ S . Moreover, \/ — A will denote the Fourier multiplier of 
symbol |£|. Observe that \J — A = D 1 . 

Let Rj be the Riesz transforms, defined via Fourier transform by Rj(f> = i 
for 3 = 1, 2. 

The divergence operator will be denoted by V- or div . 

We denote by || • \\l,p (1 < p < 00) the standard norm of the Lebesgue spaces 
L p (R n ) (n = 1, 2) and (•, •) the scalar product in I?. 

If v = («!, v 2 ) T e L 2 (R 2 ) 2 , then wc write j|vj| L 2 = (||i)i|| 2 2 + ||i' 2 i|| 2 ) 1/2 . 

If v = {v 1 ,v 2 ) T G L°°(IR 2 ) 2 , then we write ||v|| L =c = ||wi|| L = + ||w 2 ||l-- 

The standard notation H s (M. n ), or simply H s if the underlying domain is clear 

from the context, is used for the L 2 -based Sobolev spaces; their norm is written 

II -Ik- 



v 



Finally, if u = u(x, t) is a function defined in R 2 x [0, T), respectively in w 



x 



and 1 < p, q < 00, we define the mixed space-time spaces L\L V X ^ respectively L^LP, 
by the norms 

NU|XS = ( f T \H;t)\\hdt)\ resp. \\u\\ L , Ll = ( / \\u{; t)\\%dt) \ 
and L^L^p, respectively L^Ll, by the norms 

\W\\l*l*= ( \\ u (x,-)\\ P L idx) P , resp. \\u\\ l p l <, = ( \\u{x, ■)\\ p Lq dx) P . 



2. Classification of Boussinesq systems 

The Boussinesq systems can be conveniently classified according to the lineariza- 
tion at the null solutions which display their dispersive properties [7]. More pre- 
cisely, the dispersion matrix writes in Fourier variables, 
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gi(l-ec|g| 2 



■ 6(l-ee|g| 2 

\TwlfT 



6(l-ea|;| 2 ) 
l+efe|«| 2 ) 









g 2 (l-6a|g 3 | 2 ) \ 



The corresponding non zero eigenvalues are 

-(l- £a |£| 2 )(l- ec |£| 2 ) 



A±=±*|£| 



(l+ e ^| 2 )(l+66|C| 2 ) 



Recall [6] that the well-poscdncss of the linearized Boussincsq system requires 
that b > 0, d > and a < 0, c < 0, (or a = c). 

The order of A± determine the strength of the dispersion in the Boussinesq 
systems, the more dispersive one corresponding to b = d = 0, the so called KdV- 
KdV case which will be studied in details in the next section. 

By "weakly dispersive" Boussinesq systems, we mean the case where b > and 
d > 0, so that A± have order 1, or — 1. 

This situation has been studied in |15j where it is established for instance in 
the generic case with a < and c < that the Cauchy problem is well-posed in 
^([O,^];^ 1 ^ 2 )) 2 where T e is of order O(^) for any (5 < \. 

As already mentioned, existence on time intervals of order 0{\) has been estab- 
lished for weakly dispersive Boussinesq systems when the initial data are smooth 
(and with a loss of derivatives) in [27] . 

3. The KdV-KdV Boussinesq system 



(3) 



We are interested in the following dispersive system 
dtr\ + div v + ediv (?yv) + ediv Av = 



Xi,X 2 ) G 



t G 



d t v + V?y + 4 V (I V I 2 ) + eVA? 7 = 

which corresponds to the case where b = d = and a = c = i. We have scaled 
a and c to the value 1. This system might not be the best one as modeling and 
numerical purposes are concerned (in particular the cubic dispersion terms induce 
serious numerical difficulties). Nevertheless ([3]) is a mathematically interesting 
system since it can be written as a nonlocal two-dimensional version of a KdV 
-type system. Previous results concerned the one-dimensional version: 



(4) 



Vt + v x + e(?7«)x + ev xxx = 

v t + Vx + ^{v 2 )x + < n, , , = o 

Actually, as noticed in [7J , the change of variable rj 
§4$ to the following system: 



x G R, t G 



u + w, v = u — w r 



educes 



(5) 



ut + u x + eu xxx + e[2uu x + [uw) x ] = 
Wt — w x — < a-, , , + e[~2ww x + (uw) x ] = ' 

which is a system of KdV type with uncoupled (diagonal) linear part. Thus (see 
[7J) the Cauchy problem is easily seen to be locally well-posed for initial data in 
iP(R) x H S (R), s > | by the results in gD], [2T]. On the other hand, as noticed in 
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[29j Appendix A in a slightly different context, a minor modification of Bourgain's 
method as used in [23] allows to solve the Cauchy problem for ([5]) for data in 
H S (R) x H S (R) with s > — |. We refer to [H] for details. It is worth noticing that 
in [8] the question of the dependence of the existence time with respect to e is not 
considered. 

Coming back to the two-dimensional system, wc will establish that the Cauchy 
problem is locally well-posed for data in H S (R 2 ) x H S (R 2 ) 2 where s > |. 

Theorem 1. Let s > | and < e < 1 be fixed. Then for any (r] , v ) £ H S (M. 2 ) x 
H S (R 2 ) 2 with curl Vq = 0, there exist a positive time T — T(||(?7o, Vo)||_£p x (iP) 2 )> a 
space such that 

(6) Yf € C([0,T e ] ; H S (R 2 ) x H S (R 2 ) 2 ), 

and a unique solution (t],v) to ([3]) in Yfi satisfying (77, i?)| t=0 = (^Oj^o); where 
T t = Te~i. 

Moreover, for any T' G (0, T e ), there exists a neighborhood fl s of (t)o,Vo) in 
H S (M. 2 ) x H S (S. 2 ) 2 such that the flow map associated to ^ is smooth from Q, s into 

Remark 1. Our proof applies as well with minor modifications to Boussinesq sys- 
tems with b = d = and a < 0, c < 0. While this case is excluded for systems 
modeling purely gravity waves, it may happen for capillary-gravity waves when the 
surface tension parameter r is large enough. 

Remark 2. It transpires from the proof of Theorem Q] that the solution is uniformly 
bounded (with respect to e) in the corresponding spaces on a time interval [0, ^=], 

The strategy is first to diagonalize the linear part of ©. We thus reduce the 
original system to a nonlocal one. In particular, the nonlinear part involves order 
zero pseudo-differential operators (in fact Riesz transforms). We then solve the 
underlying Duhamcl integral formulation by a fixed point argument in a ball of 
a Banach space constructed from the various dispersive estimates satisfied by the 
linear part. Note that here, we are able to compute the dependence on e of the 
constants appearing in the linear estimates. This allows us to obtain an existence 
result in a time interval [0,T] depending on e, in our case T ~ e~2. 

We will thus proceed as follows. We first derive the Hamiltonian formulation of 
((3|) . We then perform the diagonalization and state the dispersive estimates which 
are to be used in the fixed point argument. The maximal function and Strichartz 
estimates seem to be new. Finally, we solve the new system (|10[) which is equivalent 
to ((3]), assuming curl Vq = 0. 



3.1. Hamiltonian structure. The system ([3]) can be rewritten on the form 

(7) d t u + A c u + ejV(u) = 0, 

where 

(1 + eA)8 Xl (1 + eA)8 X2 

(l + eA)8 Xl 

(l + eA)d X2 
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and 

Af(u)=[ \d xi {vl + vl) 

We will denote by (■,■) the scalar product on L 2 (R 2 ;M 3 ), i.e 
(u, u) = / (-qrj + viv\ + v 2 v 2 ) dx\dx 2 

JR 2 

and by J the skew adjoint matrix 





^0 


d Xl 


d X2 














\px 2 









Then, the system §5§ is equivalent to 



/(l + eA)77+ f|v 
d t u = — J (1 + eA)«! + erivi 
\ (1 + eA)v 2 + erjv 2 . 



J(gradff e )(u), 



where H e (u) is the functional given by 
1 



H e (u) 



e\Vr]\ 2 + e|Vv| 2 -if 



e?7|v| 2 )(ia;i(ia;2. 



Therefore, it follows that H e is a conserved quantity by the flow of ([3]), since 

j t H 6 (u)=H' e (u)d t u= ((gradff e )(u),9 t u) = ((gradtf e )(u), J(grad# e )(u)) = 0, 
where we used the fact that J is skew adjoint. 
3.2. Linear estimates. 

3.2.1. Diagonalization. We transform here ([3]) into an equivalent system with a 
diagonal linear part. 

First, we observe that the Fourier transform of A is given by 

/ *(l-e|£| 2 )a z(l-6|e| 2 )6 N 

A e (£i,&)=kl-e|e| 2 )£i 

W-e|£| 2 )6 

We compute the characteristic polynomial of A^i,^), 

XA e( e l! s 2 )( A ) = - A ( A2 + ( 1 -^l 2 ) 2 ^ 2 )' 

so that its eigenvalues arc given by 

A =0, Ai = i(l - e|£| 2 )|£|, and A 2 - - e|£| 2 )|£|, 
with associated eigenvectors 





( o \ 








E = 


Ss. 


E, = 




and E 2 




\ \i\J 
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Then, if we denote by P the matrix of {Eq, Ei, E 2 ) in the canonical basis, we have 
that 



P 



( 





1 








£2 


Si 




\ 


\f 




\- 















TET LET LET I and P 



( 
1 



if 



\ -1 

V 1 i«i 



Therefore, the linear part of is equivalent to <9tW + Dw = 0, where 

1 




P- 1 = 




and 



D 



'0 
i(l + eA)\/^A 



^0 -i(l + eA)V-A, 

Next, we turn to the nonlinear part of 0. Af(u) is given as function of w by 



( dx\ ((wi - w 2 )(R 2 w + Ri(wi + w 2 ))) + dx 2 (K - w 2 )(-Riw + R 2 (wi 
\d xi ({R 2 w + RxWi + Riw 2 ) 2 + {-Riw + R 2 w x + R 2 w 2 ) 2 ^ 



2 U X2 



(R 2 w + i?iu>i + R\w 2 ) + (—Riw + R 2 wi 



Then, deduce, using the identities R 2 d xi 
that 



P~W(u) = 41 I + I. 




w 2 ))) N 



-A, 



where 

/ = I(w 1 ,w 2 ) ={wi~w 2 )V^A{w 1 +w 2 ) + d Xl (w 1 -w 2 )R 1 (wi+w 2 ) 

(8) 

+ 3 X2 (wi - w 2 )R 2 (wi + w 2 ), 

and 

(9) II = II(w u w 2 ) := i\/ = A ((PiOx + w 2 )) 2 + + w 2 )f) . 

Summarizing we have that ([3]) is equivalent to 



(10) 



d t w 1 +i(l + eA) v / ^Aiv 1 + (I + II){w u w 2 ) = 
d t w 2 - i(l + eA)-/^^ + {-I + II){w u w 2 ) = 



where I(wi,w 2 ) and II(wi,w 2 ) are defined in ((Sj> and (j9|). 

Remark 3. Note that we use in our analysis that wq = 0. Indeed, the equation on 
Wo is dtWo = 0. Moreover, 



w = 



R2V1 = i?1^2 



curl v = 0. 



We observe that this condition is physically relevant. The Boussincsq systems de- 
rived from the water waves equations, where the fluid is supposed to be irrotational 
and v is an 0(e 2 ) approximation of the horizontal velocity at a certain depth which 
is a gradient. Note also that since the equation for v writes dtv = VF, the condition 
curlv = is preserved by the evolution. 



TWO-DIMENSIONAL SURFACE WAVES BOUSSINESQ SYSTEMS 



9 



(11) 



Now, to derive the smoothing effects associated to the iinear part of ([3]), it 
suffices to consider the linear system 

Ut ± L e u = 0, 
m(.,0) = u 0> 

where L e = i(I + eA) \/— A = —iLp e (D), i.e. <p e (£) = e|£| 3 — In the sequel, we 
will identify (p e (£) = (p e (\£\). Moreover, we will denote by Uf(t)uo the solution of 
CEU), i.e. 

(12) U?(t)uo= (e ±lt ^y. 
The smoothing effects are of three different types. 

3.2.2. Dispersive smoothing estimates. We will use the results in |19j (see also the 
general results in |25j ) to deduce the following local smoothing estimates. 

Let {Qa} QgZ 2 denote a family of nonoverlapping cubes of unit size such that 
M 2 = U Q a . 

Theorem 2. Let T > and e > 0. Then, it holds that 

(13) supf/ / |P - h Dlut{t)u {x)\ 2 dtdxY <c-i||«o|Uj, 

a K JQ cl Jo >c 7 

where the implicit constant does not depend on e and T . 

Proof. Without loss of generality, we can assume that Q a = Q := {x : \x\ < 1}. 
We fix a function r\ in Cq°(R) such that < rj < 1, f] = 1 on [—1,1] and suppr? C 
[0,2]. Let us define rj e {x) = e^r](e^x). Then P < _i , respectively P _i denote the 
operators defined by 

(14) P^ e _ iUo =^- 1 ( ?7e (|.|)So) and P >e _i =1-P^_j. 

Observe that in the support of 1 — r] e (\ ■ |), we have that |V<£> e (£)| = 3e|£| 2 — 1 > 0. 
Then it follows by using Theorem 4.1, formula (4.2), p 54-55 in []j5] that 

\\DlUt(t)P > e -HL= £ ( / ! q if i Mg)| 2 ^ 

(15) >e "^qx[o,t] \J\z\ >e -i 3e|£| 2 - 1 

£ e"^ll"o|U|, 

which concludes the proof of Theorem [2J □ 

3.2.3. The maximal function estimate. We will prove our maximal function esti- 
mate in the n-dimensional case. In other words, we will consider the unitary group 
Uf{t) = e ±lt ^^ D \ where p e (f) = p 6 (|f|) = e|£| 3 - |f| and £ e E n , for n > 2. Let 
{Qa}aez n denote the mesh of dyadic cubes of unit size. The main result of this 
subsection reads as follows. 

Theorem 3. With the above notation, for any s > ^ , e > and T > satisfying 
eT < 1, it holds that 

i 

(16) (V sup sup \U±(t)u {x)\ 2 Y < (1 + tS-1)\\uo\\h; 

K a&n \t\<Txe Qa > 

where the implicit constant does not depend on e and T. 
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We will only treat the case of Ul, since the case of Ut is similar. The proof of 
Theorem [3] is based on the next lemma. 

Lemma 1. For k e Z +; let t/j k £ ([2 k - 1 ,2 k+1 }) be such that < tp k < 1- Then 
for et£ (0,2], 

(17) / e^^+^Vfcdei)^ <cH k {\x\), 
where H k is decreasing and satisfies 

(18) / H k (\x\)dx<(l + \t\*-t)2 s &, 



andH k (r) < c2^ forr £ (0, 10). Also, a similar result holds for ip £ Cq°([— 10, 10]) 

3kn 

with 2~2~ replaced by c. Observe that the implicit constant does not depend on e 
and t, but may depend on the dimension n. 

The following estimate, essentially proved in Proposition 2.6 of [2D], will be 
useful. 

Proposition 1. For k £ Z+, et £ (0, 2], r £ R and Vfc as * n LemmaUi define 

r+oo 

(19) 4(i,r)= / e^ ts +sr ^Ms)ds. 

Jo 

Then 

c2 k for \r\ < 1 

(20) \h(t,r)\ < F k (r) = { c2tr"5 for 1 < \r\ < c2 2k 



cr~ N for \r\> c2 2k 



for any N £ 



Moreover, it is known (see [32] for example) that the Fourier transform of a 
radial function /(|a;|) = f(s) is still radial and is given by 



(21) f(r) = /(|e|) = r"T- / /( s )J^(r S ) S = ds, 

Jo 

where J r „ is the Bcssel function, defined by 

(22) J m (r) = <>V^r__ ^ e irs(i_ s 2jm-i d ^ for m> 



T(m + 1/2)tt2 7-i 2 

Next, we list some properties of the Bessel functions (see [32], [22], [16] and the 
references therein). 

Lemma 2. If ZioZds that 

(23) J ro (r) - 0(r m ), 

r— ¥0 

iV JV 

(24) J m (r) - e" lr Va mj T-^' + 5) + e ' r V5 m ,r"«+5' + 0(r-( w +l>), 

r— >+oo L — ' * — J 
3=0 3=0 

for any N £ Z+, and 

(25) r-^J^(r) = c n U(e lr h(r)) , 
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where h is a smooth function satisfying 

(26) \d^h(r)\<c k (l + r)-^ 1 - k , 

for any k £ Z + . 

Proof of Lemma [0 By using (|2~Tj) , it follows that 
(27) 

X k (t,r):= / e^'^ +x ^M\^\)dC = r-^ / e ui -""- s ^ k {s)J^2{rs) S ^ ds, 



where r =\x\ £ [0, +oo) and et £ [0, 2]. 

When 0<|r-£|<lor0<r< 1, it follows from (gSJ) (or (Hn])-®) that 

r+oo 

(28) \X k (t,r)\< i> k { s ) s n - 1 ds = c2 nk . 

Jo 

In the case r > 1 and |r — t\ > 1, we substitute J n-2 by the right-hand side 
of (|24[) in (|27p and evaluate successively each term of the sum and the remainder. 
Here we consider only the most difficult case, when a m ,j = 0. Then, the j th term 
has the form 

l kjj {t,r) :=r-^ / e l(t{ts3 - s)+sr) 4> k (s){rs)- {l+ ^ a* ds, 



r+oo 

\T k ^{t,r)\<r-^{2 k r)-^+^2 kl i / e^ tes3+s ^Ms)di 



so that 



where tp k is another function satisfying ip k £ C 3O ([2 fc 1 ,2 k+1 ]) and < tp k < 1. 
Therefore, we deduce from Proposition [T] that 



(29) 



IT (t ^l</ 2 fc( t- J V-(t-Hj)|r-t|-5 for \r-t\ £ [l,c2 2k ] 

~ | 2fe ( f _ J -_i) r _ (f _i +j )| r _ t |_ m for | r _ f | >c2 2fc 5 



for any m £ Z + . Next, we fix N = N(n) > ^^—^ and bound the remainder 

/>+oo 

H k (t,r) :=r-^ / e titip ' {s)+sr ^ k (s)(rsy iN+ ^s^ds 
Jo 

as follows 

(30) \n k {t,r)\ < 2 k ^- N ~^r"^ +N+ ^ < r~ m , 

with m > n. 

Therefore, if we define 

( 2 kn for \x\ < 1 or ||a;| - t\ < 1 

Wfc(N) = < E^o 2fc(5_i) M~ (f ~' +i) N for 1 < IN -*l < c22fe 

[ 2 k ^-^\\x\-t\~ m for ||x| -i| > c2 2fe , 

and 



ff fe (|z|) =W fc (|a;|) + 
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it follows from (|27l) - (f3"u]) that \lk{t, \x\)\ < cHk(\x\) and a simple computation leads 
to 

/ H k {\x\)dx < c2 kn + c2^ [ r%-*\r-t\-*dr 

JK" Jl<\r-t\<c2 2k 
< C 2^(l + 

which concludes the proof of Lemma [TJ □ 

Finally, at this point, the proof of Theorem [3] follows closely the argument of 
Kenig, Ponce and Vega in the case of the Schrodinger equation in Theorem 3.2 of 
l22l. Therefore, we will omit it. 



3.2.4. Strichartz estimates. Strichartz estimates for unitary groups of the form 
e *t4>(D) m jgm^ n > 2, were derived in the case where <j) is an elliptic polynomial 
by Kenig, Ponce and Vega [19j and in the case where is a general polynomial in 
M 2 by Ben-Artzi, Koch and Saut [3]. When the phase function is a radial (non- 
homogencous) function and its derivative does not vanish, some techniques were 
recently developed by Cho and Ozawa [H] and Guo, Peng and Wang [TB]. 

In the sequel, we will use the techniques developed in [TB], based on the ones 
used in |19j and on the representation of the Fourier transform of a radial function 
in terms of the Bessel function (see formula (|2"Tjl ). to prove Strichartz estimates 
associated to the unitary groups Uf- defined in (fT2|) . However, in our case, we do 
not need to perform a dyadic decomposition in frequencies. 

Theorem 4. Let < e < 1, T > and < a < |. Then, it holds that 

(31) \\D2U?uo\\l<°ls> ^^-IMI^, 

for all uq G L 2 (R 2 ), where the implicit constant is independent of e and T , q a is 
the root of the polynomial 

3q 2 - 2(7- 2a)q+ 12 = 0, 
satisfying q a > 2 and K a = ^ + ^ — . 

Remark 4. When a = 0, then q = 7+ f^ = |+ and K = \ - = f +. On the 
other hand, we have that lim Q _ i .i q a = 2 and lim^^^^ K a = | 

For sake of simplicity, we will fix U e = in the rest of this subsection. First, 
we derive the following decay estimate for the solution to the linear problem 

Proposition 2. Let < e < 1 and < j3 < 1. Then, it holds that 

(32) \\DPu e (t)uo\\ Ls > <k p>£ (t)\\u \\ L i, 
for all f 6 I, where kp, e is given by 

if t < e£ 

if t>6ei ' 
and 6 is any positive constant independent of e. 



(33) kp, e (jt) 




The proof of Proposition [2] is based on formula ([21]) , Lemma [2] and Van der 
Corput's lemma: 
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Lemma 3. Suppose that f is a real valued C 2 '-function defined in [a,b] such that 
|/" (01 > 1 for any £ G [a, b]. Then 

where the implicit constant does not depend on a and b. 

Proof of Proposition [H First observe that 

(34) WDfofluoh- < ||(l • fe-^ILJMU 1 - 

On the other hand, formulas ([21]) and (|25[) imply that 

(35) 

/"+00 

(|e|^e- ii ^) V (x)= / ^e^'-^JoMatfa 
Jo 

/■+oo r-\-oo 

= / / e ,f(£s3 - s) e irs Mn)s(is+ / s^'^-^e-^/^Jsds 
Jo Jo 

= Ip{r)+II p {r) 

where r = \x\. 

For sake of simplicity, we will assume that t > and only deal with Hp since 
can be handled by similar techniques. We change variables u — (ei)ss and deduce 
that 

(36) II (r) = (et)- 2 -^A p (-^- T ), 

(et) 3 

where 

r+oo 

A p (r) = / e ifr ^h(ru)u 0+1 du, 
Jo 

the phase function f r (u) is given by 

/ r (ii) = u 3 — (r + a)u, and a = t^e~^. 

Therefore 

(37) sup \A p (r)\ < max{l,a7+3} = max{l,t^ + h~^~^}, V/3e[0,l], 

r>0 

would imply formula (|32[) . 

To prove (|37|) , we introduce the smooth real- values functions (ipi , 1P2) G x C 00 
such that < Vii V>2 < 1, i>i(u) + "02 ( u ) = 1, 

supp-01 C {u : |3u 2 — (r + a)\ < "t 

and 

^ 2 =0 in {u : |3u 2 - (r + a)| < — — }. 

O 

It follows that 

(38) \Ap(r)\<\A}(r)\ + \A}(r)\, 
where 

r+oo 

AUr) = / e lfAu) h{ru)u p+1 i) 3 {u)du, j G {1,2}. 
Jo 

First, we deal with ^4 2 . Observe that we can restrict the oscillatory integral 
in the range u G [l,+oo], since otherwise when u G [0,1], the estimate is trivial. 
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Moreover, we deduce from the triangle inequality that in the support of ifo, the 
derivative of the phase function satisfies \f' r {u)\ = |3u 2 — (r + a)| > \{u 2 + (r + a)). 
Then, we obtain integrating by parts that 

J i ij r {u) au 



du V f' r (u) J 

Therefore, it follows from (|26p that 
(39) 

f + °° / riiP +1 iif 3 + 2 \ 

\Al{r)\ < / f r + r )du < 1. 

A v (?i 2 + (r + a))(l + ru)2 (v? + (r + a)) 2 (l + ru) 5 ' 

Note that the implicit constant does not depend on r, e or t. 

Next we turn to Al. In the support of ip%, we have u ~ (r + a) 2, so that 

I /"(u) I = 6u > (r + a) 3. Thus, it follows from Van der Corput's lemma that 

^ (r + a) i V aw 

(40) (r + a)hW) .. £ + 

< ^ j- < max{l,a 2 + 4 }. 

(r + a) 4 (1 + r(r + a) 2 ) 2 

Finally, we deduce formula (f3"T|) combining (f55|) . which concludes the proof 
of Proposition [2] □ 

We are now in position to give a proof of Theorem 01 

Proof of Theorem^ Fix < a < i, j3 = 2a G [0, 1), k = n a , q = q a and q' its 
conjugate exponent, i.e. - + ^7 = 1. We first observe by using a P. Tomas' duality 
argument (see for example [2S]) that estimate (|3T|) is equivalent to 

(41) 11/ D0U e (t-t')g(;t')dt'\\ LqLoo <e- 2 «\ 



for all 3 G L q ' (M.; L 1 ^ 2 )). 

Next we prove estimate (jUJ. It follows from Minkowski's inequality and estimate 
(|52|> that 

/•+OO 

(42) 11/ Df[/ e (i -*')<?(•, *||ff(-)IU;(*) :=^,e(*), 



where fc^ i£ is defined in (j3"3")l . We will denote <p(t) = ||g(-, i)||ii ■ Then we divide 
the kernel fc^ j£ in two parts, kp^ e = A;2 e + fcej° e , where 

*g, e (0 = W- 2 -^x {meei} and ^(i) = rH r i X{|i|yii}I 

so that 

(43) j^(t) = 4Jt) + J™M> 

where J° e , respectively J^° e , is the convolution operator associated to the kernel 
hp e , respectively fc^° e . 
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To estimate J2 e , we observe that 



/ k° 0e (t)dt = 2 / (et) _a t-dt = 



2 

2+£ , 1+0 1-3 

ce 2 w 3 



since < ft < 1. Therefore, it follows from Theorem 2 in Chapter III of [21] that 
(44) \J^(t)\<e^6^M^t), 

where M. denotes the Hardy-Littlewood maximal function. On the other hand, 
Young's theorem implies that 



(45) \J% e (t)\ < ( [ t t-idb) |M|jy =c 



= «-< 1+ *-*>**-*|M| I y 1 



'|t|>9e' 

J 

2 2q ^ 2 



since g > 2. Observe that 1 + f — ^ > Thus, we deduce gathering (|33" ]) -([IS |) 
that 

(46) |Jp, e (t)| <e- {1+! *-^\0°-?\\<p\\ L «> +6^Mip{t)). 

Now, we choose 6 = 6{t) to minimize the term on the right-hand side of (|46|) , which 
is to say 

6{t)i-^ = \\<p\\ Lq ,My{t)-\ 
This implies together with (|4l)]) that 

(47) |J^,e(t)| < e - (1+ ^ ) |M|^,.M^) 1 -'', 
where 7 = 5^2/^-6 • Then it follows that 

(48) ll^elU^e-^-^IK^HMvll^- 
Moreover, observe that 

(l-7)? = 8 / 3g 2 - 2(7-/3)^+12 = 0. 

Finally, we conclude from the fact the maximal function is continuous in L q that 

(49) \\JeAi*<<r^-&\\<p\\ L <, 

which concludes the proof of Theorem 01 □ 

3.3. The nonlinear Cauchy problem. The objective of this subsection is to 
prove the following result. 

Theorem 5. Let s > | and < e < 1 be fixed. Then for any (w^w®) € H S (R 2 ) x 
H S (M. 2 ), there exist a positive time T = T(\\(wi,W2)\\h s xH s ), a space such 
that 

(50) X s Tr ^ C([0,TJ ; H S {R 2 ) x H S (R 2 )), 

and a unique solution (w\,W2) to (|10[) in Xj> satisfying (w\,W2)\ t=Q — (w^w®), 
where T c = Te~? . 

Moreover, for any T' € (0,T e ), there exists a neighborhood f2 s of (iuJ, u^) in 
H S (R 2 ) x H S (R 2 ) such that the flow map associated to (|10[) is smooth from fl s into 
X?pi • 

First, we list some well-known properties of the Riesz transforms. 
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Proposition 3. It holds that 

(51) \\Rjf\W < C\\f\\ L r, 
for all 1 < p < oo and j € {1, 2} ; and 

(52) Dl = Rxd Xl +R 2 d X2 . 

The following Leibniz' rule for fractional derivative, derived by Kenig, Ponce and 
Vega in Theorem A. 12 of [21], will also be needed. 

Lemma 4. Let < 7 < 1 and 1 < p < 00. Then 

(53) \\Dl(fg)-fDlg-gDlf\\ Lp <\\g\\ L ~\\Dj<f\\ LP , 
for all /, g : R" C. 

Proof of Theorem [5l We will treat only the most difficult case when | < s < 2 and 
we define 7 = s — 1 £ (|, 1). The integral system associated to (fTu| with initial 
data (wi, w®) can be written as 

(54) \ Wl =J+(wi ' W2) :=C/ e + (*K + e/o U+(t-t')(I + II)( Wl ,w 2 )(t')dt' 
\w a =F-(w 1 ,W2):=U-{t)w% + ef*U-(t-t')(-I + II)(w 1 ,w 2 )(t')dt 

where the nonlinearities / and II appearing on the right-hand side of (|54|) are 
denned in (0) and ©. 

Next, for < T < e~ 2 , we fix a = 1 — 7 G (0, ^), g Q and k q defined in Theorem 
S] and consider the following semi-norms 

A^/)= sup ||/(t)||*., 

0<t<T 

a^(/) = ^ + E ii^/n L i +ioc + el+ ii^/n i i +L . + el+ E E 11^/11^ 

|^|<i T t x ;=i,2|^|<i T 

+ ^ + E Px^/ll A +T , 
1=1,2 Lt L * 

\Uf) = e K " E 11^/11^^+^ E E \\Dl%Rif\\Lt°>L~, 

1/9 1 = 1 !=1,2|/3|=1 

= e* E ™PjP >e -i€D x + V\\ L > m 



+ E E sup HP^i^D^/l 

J=l,2 |/3| = 1 q£Z2 

4(/) = c*( E \\f\\l^ aXl0iT] ) h +^ E ( E \\ R if\\L~ aXl0tT] 



d \\™U if 

(=1,2 q£2 



where {Q Q } QgZ 2 denotes a family of nonoverlapping cubes of unit size such that 
R 2 = U Q , e z2<3a ! as in the precedent subsection. Wc also define the Banach space 
X| by 

Xl = {(w x , w 2 ) € C([0,T];H s (R 2 ) x H S (R 2 )) \ \\(w 1 ,w 2 )\\ X } <°°}- 

where 

4 4 
II (wi,w 2 )\\ x2t = E A J T K) + E A J(^)- 
i=i i=i 
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We will show that, for adequate T and a, the map (J 7+ , J~~) is a contraction in a 
closed ball X^(a) in X|> of radius a > and centered at the origin. 

Using the integral equation (JMJ), Minkowski's integral inequality, the linear es- 
timates (fT3)) . (fT6| and ([3"Tj) with a = and a = 1 — 7, and the fact that the Ricsz 
transforms are unitary operators in H s , we deduce that 

4 4 
K ( («* . w 2 )) + 51 4 ( ^" K , W 2 ) ) 

3=1 3=1 

(55) T 

< HK,^2)llH-xff. +e / (||/(«; 1 ,t«a)(t)||a. + ||J/(t£;i J ifl3)(*))||H S )dt. 



According to formulas (JSJ and ©, the nonlincarities /(u^u^) and II(wi,W2) 
contain terms of the three following types: WjD x w k , d Xl WjRiWk and D x {RiWjR m Wk) 
for j, fc, Z, m = 1, 2. For sake of simplicity, we only will present the computations 
for a nonlinear term of the first kind, since the other ones can be handled by similar 
arguments. In other words, we have to estimate the term e \\wjD x w k \\H^dt as 
a function of the norms defined in (|3.3p . By the definition of the Sobolev space 
H S (R 2 ), we have that 

(56) e[ \\w 3 Dlw k \\ Hi dt<e( \\w 3 Dlw k \\ Ll dt + e [ \\D s x { Wj Dlw k ) \\ Ll db. 
Jo Jo Jo 

First, we use Holder's inequality and the Sobolev embedding H S (M. 2 ) i°°(M 2 ) 
to estimate the first term on the right-hand side of (|56")l as 



(57) 



: / \\wjDlw k \\ L 2dt < eTWwjW^L^WDlwkWL^^ 
Jo 

< eTXip(wj)X^p(wk)- 



To treat the second term appearing on the right-hand side of (|56p . we observe 
from {55) that D% = D2R\d Xl + D%R 2 d X2 . Therefore, it follows from flSTJ that 



/ \\D x ( Wj D l xWk )\\ Ll dt 
Jo 

<^E / \\ D Z( 9 *M D l™k)hldt + eJ2 I \\Dl{w 3 d Xl Dlw k )\\ Ll dt. 

1=1.2 1=1,2 



6 

(58) 



Then, we deduce from the fractional Leibniz' rule (|53|) that 
(59) 

T 

WDHd^Dlw^dt < e^T^e^ + \\d Xl w 3 \\ r r + r jD 1 x + -'w k \\ L¥Ll 



+ er-Tr+ e y + \\Dlw k \\^ +L Jd Xl DZw j \\ L¥Ll 
< e^-T^ + (X^(w 3 )X 1 T (w k ) + \^(w k )X 1 T {w 3 )). 
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On the other hand, by using again Leibniz' rule and Holder's inequality, we 
observe that 

e f \\Dl{w 3 d Xl D l x w k )\\ Ll dt 



<e \\D2w j \\ Ili \\d asi Diw k \\ Ls ,dt + e \\w J d Xl D 1 +-'w k \\ Ll dt 



(60) 



We deal with the first term on the right-hand side of ([BOf. We deduce by using 
Holder's inequality in time that 



' [ \\ D 2i»jh*\\dxiDlwk\\L«>dt 
Jo 



(61) < e^T< |p> i |U-L S c"-||»x l I>>fc|U S - i «. 

<J-">'T%\*(w j )\Z(w h ), 

where q' a is the conjugate exponent to q a . 

Finally, to estimate the second term on the right-hand side of (|S(J]) . we observe 
that 

(62) 

< e 



\w j d xl Dl + ' r w k \\ L 2dt 



On the one hand, we use the Sobolev embedding H S (M. 2 ) c -s> L°°(IR 2 ), since s > 1, 
to obtain that 

J 

(63) <^\\ Wj \\ L ^ Hi e-h\\Dl^ Wk \\ L o 



< eiT\^(wj)\^(w k ). 
On the other hand, we deduce from Holder's inequality that 
(64) 

e / \\w 3 Pid Xl Dl^w k \\ Ll dt 
Jo 

^ eT HYl I \ w i p ^ -ld xl Dl^w k \ 2 dxdt 

\1^2«x[o,t] > £ 

< eTU V ||w,|||oc V sup IIP _id x ,D 1 x + ^w k \\ L 2 

< e f Ti\^( Wj )\^(wk)- 
Thus, it follows gathering ([5"S ]) - (j62]) that 



(65) 



\(j r+ {w 1 ,w 2 ),J 7 (w 1 ,w 2 ))\\ X 2 < c\\(w°,W2)\\h3xH* 



c(e?-T?+ + e^TZ + £ 3T + e«T*) w 2 )| |2 
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and similarly, 

||(j r+ (wi,U;2),-7 r ~(wi,U'2)) - ( J r+ ( Wl , W2 ) , ( Wi , 7l7 2 )) 1 1 ^ 

(66) < c (e?~T7+ + e 1 - K "T^ +e^T + e^T^) 

x (|| (wi, w 2 )\\xl + II (*i, w 2 )\\xl)\\{wi,w 2 ) - (wi,w 2 )\\xl- 

Observe that 

c(e?~T?+ + e 1 ^ Ra T 1 t + e ?T + e?T5) < 1, 

since T < e - '. Therefore, we deduce by choosing 

a = 2c||« W 0)|| H2xH 2 and T ~ e"? 

that (J-" + , J 7- ) defines a contraction in the space X^(a), which concludes the proof 
of Proposition [5] by the Picard fixed point theorem. □ 

4. Other strongly dispersive Boussinesq systems 

Other two-dimensional Boussinesq systems can be classified according to the 
order of the eigenvalues of the dispersion matrix (see Introduction) . As we already 
mentioned, when b, d > 0, the well-posedness on time scales of order 1/y/e has been 
established in [15] for initial data in low order Sobolev spaces and on time interval 
of order 1/e in [27] for smooth initial data (with loss of derivatives). We thus 
restrict ourselves to the case where at least one of the coefficients b or d vanishes, 
emphasizing the strongest dispersive case corresponding to eigenvalues of order 2 
(the so-called "Schrodinger" type systems) where a < and c < 0. 

The local well poscdness for some of those systems has been established in the 
two-dimensional case by Cung The Anh [T3], but without looking for the depen- 
dence of the existence time with respect to e. Note however that the proof given in 
[T5] for the case where d = (Theorem 3.5) does not seem to be correct. 

We will complete the analysis, in particular by checking the e dependence. It 
turns out that when d > 0, one can establish the local well-posedness of Boussinesq 
"Schrodinger" type systems on time intervals of order l/y 7 ^ by elementary energy 
methods, in relatively low order Sobolev spaces. 

We first state a useful lemma. 

Lemma 5. Let s > 3/2. Then there exists C > such that for any v € H S+1 (M. 2 ) 
and i] e iJ s (M 2 ), 

(67) |(A'V-(» 7 v),A^)|<q|v|| H . + iN|i r .. 

The proof of this lemma is based on the following commutator estimate obtained 
by Ponce in Lemma 2.3 of |28j . 

Lemma 6. If s > 1 and 1 < p < oo, then 

(68) ||[A S ;/]<?|U P < CdlV/II^IIA 8 "^!!^ + \\A s f\\ L , 3 \\g\\ L **), 
where pi, p4 € (l,+oo] such that 

1 1 1 1 _ 1 

Pi P2 P3 Pi P 
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Proof of Lemma Fix s > |. First observe that 

(69) (A*V • fav), A s 7?) - (A'^^i), A s ry) + (A s c^ 2 (r^ 2 ), A S V ). 
Moreover, we can rewrite the terms appearing on the right-hand side of (|69|) as 

(70) (A S 8 X] (r)Vj), A S V )= (A s ( v d X] v 3 ),A s r,)+ ([A s , v^r,, A° v ) + :r,YiL ,,. A s v ). 

We use Lemma g] and the Sobolev embedding H a (R 2 ) <-» L°°(R 2 ), for a > 1, to 
estimate the first term on the right-hand side of (|70|) as 

(71) \(A s ( v d X] v 3 ),A s v )\<C\\v\\ H s + i\\ V \\ 2 Hs . 

On the other hand, it follows from Lemma [5] with p\ = oo, pi = 2, p^ = p± = 4, 
and the Sobolev embedding Hi(R 2 ) ^ L 4 (R 2 ) that 

|([A s ,^]^r,,A s r,)| 

(72) < C(\\Vv j \\ L ~\\A'- 1 d Xt ri\\ L 2 + IIA^II^II^^II^HA^II^ 
<C||v||^ +1 h||^. 

Finally, we have integrating by parts that 

( Vj A s d Xj r],,A s V ) = —(d^vjA'mA'v), 

which implies by Holder's inequality and the Sobolev embedding H a (M. 2 ) L°°(M. 2 ), 
for a > 1, that 

(73) \(v J A'a XjV ,A'ri)\<C\\v\\ H . + 4 V \&.. 

Therefore, estimate (|67| is deduced gathering (|69| -([73 |) . □ 

4.1. Schrodinger type Boussinesq systems when d > 0. As already men- 
tioned, attention is focused here in systems where the eigenvalues A±(£) are of 
order 2. Because of the constraint on a, b, c, d the only possibilities are a < 0, c < 0, 
b = 0, d > or a < 0, c < 0, b > 0, <i = 0. We will study here the first case and 
thus consider the systems 



f 77t + V • v + e[V • (r/v) + aV • Av] = 
1 j \ v t + V77 + e[iV|v| 2 + cVA77-dAv t ] =0, 

where a < 0, c < 0, d > 0. 

Those systems are called in [7] "Bcnjamin-Ono type systems" since the dispersion 
matrix has, in the one-dimensional case and when e = 1, eigenvalues equal to 
±((f )ik\k\ + r{k)) where r{k) is 0(-r|r) as k -> oo. 

In the two-dimensional case, the dispersion matrix in Fourier variables reads 



o / 

The corresponding nonzero eigenvalues are 



/ &(i 

gi(i- ec |g| 2 ) n 

l+«2|«l 2 U 

g 2 (l- ee |g| 2 ) n 



When surface tension is large enough so that r > 1/3, we have the admissible system a = b = 
d = 0, c = t — 1/3 but we will note consider it here. 
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(l-ea|g| 2 )(l-ec|e| 2 ) y 
l + ed\tf J ■ 

For e fixed, A± are equivalent as |£| — > oo to ±e (^) 2 |£| 2 , which explains the 
terminology "Schrddinger type" Boussinesq systems. 

After diagonalization of the linear part, (|74[) can thus be written as a (nonlinearly 
coupled) system of two Schrddinger type nonlocal equations involving derivatives 
of the unknowns, and we might think of solving the Cauchy problem by applying 
the general results of [21] ■ We will refrain to do that here because an elementary 
energy method on the original formulation (|74|l , which is the extension to the two- 
dimensional case of the corresponding one-dimensional result in [7] applies when 
d > 0, and when d = for curlfree velocities. 

Theorem 6. Let s > 3/2 and < e < 1 be given. 

(i) Assume that a = c. Then, for every (?7o,v ) G H S (R 2 ) x H S+1 (R 2 ) 2 , there 
exist T = T(| 1 770 1 \h s , \\'Vo\\h'+ 1 ) and a unique solution 

(r?, v) G C([0,T e ];H s (M. 2 )) x C([0, TJ; H S+1 (R 2 ) 2 ), 

where T e = TtT 1 ! 2 , of (|74|) with initial data (770, Vo), which is uniformly bounded 
on [O.TJ. 

fn) Assume that a ^ c. Then, for every (770, v ) G i7 s+1 (K 2 ) x ff s+2 (R 2 ) 2 , Tj/iere 
ea;jsr; T = T(||77o||ijs+i, ||vo||ijs+2) cmd a unique solution 

(r), v) G <7([0,T e ];iP +1 (K 2 )) x C([0, T e ]; ff s+2 (R 2 ) 2 ), 

where T e = Te -1 / 2 0/ (|74[) wii/i initial data (770, vo) which is uniformly bounded on 
[0,T e ]. 

Proof, (i) Uniqueness in the class C([0, T]; H 3 / 2 (R 2 )) x C([0, T]; ff 5 / 2 (R 2 )) 2 . 

In what follows we will use freely the embedding H 2 (R 2 ) «-» L 4 (M 2 ) and thus 
i?l(]R 2 ) W 1 ' 4 (M 2 ) (see 0]). 

Plainly it suffices to consider e = 1. Let (771, vi), (772 , V2) two solutions in the 
class C([0,T];iJ 3 / 2 (]R 2 )) x C([0, T]; iJ 5 / 2 (R 2 )) 2 and let TV = 771 - 772, V = Vi - v 2 . 
We have therefore 




(75) 



N t + V • V + e[V • (iVvi + 772 V) + aV • AV] = 

V t + VA + e[|V(V • (vi + v 2 )) + cVAA - dAV t ] = 0, 



We take the L 2 scalar product of the first equation by \c\N, of the second one 
by \a\V and add the resulting equations. The contributions of the cubic dispersive 
terms cancel out and we find 



ld_ 

2dt 



[ {\c\\N\ 2 + \a\\V\ 2 + d|VV| 2 ) = f (fliVV-V + cViV-V) 

it 2 it 2 

/ V • (JVvi + 772 V)N - \ f V(V ■ ( Vl + v 2 )) • V. 
Jr 2 1 Jr 2 



(76) 

One has successively 



One term has to be interpreted as a 

_ H-l/2 duality. 
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(77) 



\a\NW ■ V+ \c\WN-V) 



[ (\a\-\c\)NV-V 

Jr 2 



<c||iv|U 2 ||vv|| i2 . 



Then 



/ V ■ (Nvi + r] 2 V)N = [ (iV 2 V-vi+ 77 2 iVV-V + (ViV-v 1 )iV+(Vr/ 2 -V)iV. 

Jr 2 Jr 2 

Observing that J R2 (VAT ■ Vi)N = — I L 2 iV 2 V • vi, the contribution of this term 
is upper bounded by 

C(||v ■ vi|| L =||iV||| 2 + \\m\\ L -\\N\\ L 2\\vv\\ L 2 + ||V|| LJ ||v^|| L «i||iv|| L 2) 



<C(||iV||| 2 + ||VV|| 2 2 ). 



Finally, 



/ V(V-(vi+v 2 ))-V 

Jr 2 



<G(||vi|| W i.« + ||v a || W r 1 , 

The result follows from (17611 and Gronwall's lemma. 



|V||^. 



(ii) Existence. We will focus only on the derivation of the energy estimates check- 
ing the dependence of the existence interval on e. The complete proof would result 
from a standard compactness argument implemented on a regularized version of 
the system (for instance by truncating the high frequencies in the differential op- 
erators). The strong continuity in time and the continuity of the flow would result 
from a standard application of the Bona-Smith trick. 

Energy estimates in the case a = c. We apply the operator A s to both equations 
of (f74|) . multiply the first one by |a|A s ?7, take the scalar product of the second one 
by |a|A s V, integrate and sum. All the linear terms cancel out and it remains to 
estimate 



h = e / (A S V • (7?v)A s 77 
Jr 2 

Lemma [5] implies that 



1 



A s V|v| 2 -A'v) =e{I 2 +h). 



\h\<C\\w\\H^\H\ 2 H s. 
Finally, by the Kato- Ponce commutator lemma f|18j). 

N<C||v|| HS+1 ||v|| 2 ffs . 

Combining those estimates leads to the differential inequality 

(78) j t [H\ 2 H s + \\v\\%. + e\\v\\% s+1 ] < CeOlvll^lMI 2 ,, + ||v|| 2 H3 ||v|| ffs+1 ). 

If we define Y(t) = \\rj\fjjs + 1| v|||j- s +e||v|| 2 c/3+ i, then inequality ((78)) can be rewritten 
as 

Y'(t) < Ce 1/2 Y(t) 3/2 , 
which in turn implies an a priori bound on a time interval [0, where T = 
T(\\vo\\h>,\Wo\\h*+i)- 
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Energy estimates in the case a 7^ c. We apply the operator A s to both equations of 
(|74|) . multiply the first one by A s rj + ceAA 3 rj, take the scalar product of the second 
one by A s v + aeAA s v, integrate and sum to obtain that 

(79) 

~ [Nllr. - ce\\Vr,\\%. + ||v|| 2 ffs + (d- a)e\\Vv\\%. arf € 2 ||A v ||^] 
V • A s vA s n + ce I V • A s vAA s ?7 + / VA s 7y • A s v 



ae J VA s rj ■ AA s v + ae J V • AA s vA s ^ + ace 2 J V • AA s vAA s ?] 
ce y VAA s jj • A s v + ace 2 J VAA s ij • AA s v + e J A S V • {^)A s i] 

J A S V • (r/v)AA s r ? + I /" A s V(|v| 2 ) ■ A s v + ^ J A s V(|v| 2 ) • AA S 



e 2 c 



12 

Now, observe integrating by parts that J\ + J 3 = J 2 + J7 = J4 + J5 = Jg + J« = 0. 
Moreover, we deduce from Lemma [5] that 

(80) |J 9 |<e||v||^ + i|M|l, s , 
and 

(81) |Jio|<e 2 ||v|U s+2 ||7?|| 2 ff , +1 . 

On the other hand, it follows integrating by parts and using Lemma |6] that 

(82) iJiil^HvH^ilvll 2 ^, 
and 

(83) I^i2|<e a ||v|| ir .+>||v||^ + i. 
Finally, let define Y by 

Y(t) = \\ v \\* s . cllVr/H^ + M 2 H , +(d- a) e ||Vv||^ - ade 2 \\ Av|&.. 
Then it follows gathering ([79]) - (|83) that 

Y'(t) < Ce l l 2 Y(i?' 2 , 
which in turn implies an a uniform in e priori bound on a time interval [0, -^=\. □ 

4.2. Schrodinger type Boussinesq systems when d = with curl free ve- 
locity field. We thus consider the case a < 0, c < 0, b > 0, d = 0. 



(84) 



rjt + V • v + e[V • (77V) + aV • Av - bAi] t ] = 
v t + V?y + e[±V|v| 2 + cVAr?] = 0. 

We will here restrict the velocity v to irrotational motions, a situation which, as 
we already indicated is relevant since the variable v in the Boussinesq systems is 
curlfrcc up to a 0(e 2 ) term. Note also that the curlfree condition is preserved by 
the evolution of When v is curlfree, the term ^V|v| 2 in the second equation of 
(|84[) writes as two transport equations namely (v-Vfi, v-Vi>2) T where v = [v\, V2) T ■ 
This will permit to perform energy estimates on v. 
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Theorem 7. Let s > 2 and < e < 1 be given. 

(i) Assume that a = c. Then, for every (r?o,vo) G H S+1 (R 2 ) x H S (R 2 ) 2 with 
curl Vo = there exist T = T(||?7o||_H" s + 1 j ll v o||if s ) an d o- unique solution 

(77, v) G C([0,T e };H s+1 (R 2 )) x C([0,T e }; H s (R 2 ) 2 ), 

where T t = Te -1 ' 2 , o/ (|74p u>ii/i initial data (j7o,Vq), which is uniformly bounded 
on [0,TJ. 

(m) ^sswrne t/iat a ^ c. TTien, /or every (rj ,v ) £ # S+2 (K 2 ) x iJ s+1 (R 2 ) 2 with 
curl Vo = i/iere exist T = T{\\rjo\\jjs+2, ||vo||#s+i) and a unique solution 

(r?,v) G C([0,T,};H s+2 (R 2 )) x C([0, TJ; H S+1 (R 2 ) 2 ), 

where T t = TtT 1 ! 2 of (|74[) with initial data (770, vo) which is uniformly bounded on 
[0,TJ. 

The proof of Theorem [7] is very similar to that of Theorem [5] and we only sketch 
it. We will use the following estimate. 

Lemma 7. Lei s > 2 and v G i2" s (R 2 ) smc/i i/jaf curl v = 0. Then 



A s V|v| 2 • A s v 



<C\\v\\h 



Proof. Let v = (v\,V2) T - Then d Xl V2 = d X2 v\ and 

^V|v| 2 = (vid Xl Vt +v 2 d X2 v 1 ,v 1 d Xl v 2 + v 2 d X2 v 2 ) T . 

The estimate then follows immediately from integration by parts, the Kato-Poncc 
commutator lemma and the embedding H S ^ 1 (R 2 ) L°°(K 2 ). □ 

Proof of Theorem^ (i) Uniqueness. Again we consider e = 1. Let (?/i,vi), (772^2) 
two solutions and let N = 771 — 772, V = Vi — v 2 . We have therefore 

N t + V ■ V + e[V • (Nvi + 7 ?2 V) + aV • AV - bAN t ] = 
V t +VN + e[±V(V • ( Vl + v 2 )) + cVAN] = 0. 

We multiply the first equation by —cN, take scalar product of the second one by 
— aV, integrate and sum. One obtains after integrations by parts 

r |c|AT 2 + H|V| 2 + |c|6|VAf] 



(85) 



2dt 

< CdlvilU-, |MU~, HVvxIlioo, ||Vv 2 || L »)[||iV||i 2 + \\VN\\ 2 L2 + ||V|| 2 2 ], 
and the result follows from Gronwall's lemma. 

(ii) Existence. We just indicate how to obtain the energy estimates in the case 
a = c. We apply the operator A s to both equations of (|T4"]) . multiply the first one 
by A s 77, take the scalar product of the second one by A s v, integrate and sum to get 
using Kato- Ponce Lemma and Lemma [7] 

|[!|A s 77|| 2 2 + ||A s v||| 2 + e &||VA s 77|| 2 2 ] 
( 86 ) <Ce [ A S V- (77v)A s ?7 + A s V|v| 2 • A s v 

<C£(||Vv||^||A s 77|| 2 2 + ||v||^), 



TWO-DIMENSIONAL SURFACE WAVES BOUSSINESQ SYSTEMS 



25 



from which one obtain the expected a priori uniform estimate on a time interval 
[0,T/y/e\- 

When a^c one has to modify the proof as in Theorem [5] □ 

5. Final comments 

To keep this paper short and to maintain a certain unity, we have refrain here 
to consider the other Boussinesq systems, which have eigenvalues A± of order 1 or 
less. For instance the cases a = c = b = 0, d = 1/3 or a = c = d = 0, 6=1/3 
are version of Boussinesq original system. The former case has been studied in the 
one-dimensional case by Schonbek [30] and Amick [2] who obtained global well- 
poscdness by viewing the system as a (dispersive) perturbation of the Saint- Venant 
(shallow-water) hyperbolic system. 

In both cases, no results seem to be available in two-dimensions. 
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